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PREFACE 


This  report  is  part  of  RAND's  continuing  interest  in  the 
detection  of  high-altitude  nuclear  explosions.  In  particular,  it 
is  concer:>ed  with  the  low-frequency  hydromagnetic  waves  which  are 
generated  by  such  e]q>lo8ion8,  and  which  are  detected  by  ground- 
current  measurements  at  world -wide  locations.  In  this  report,  we 
direct  our  attention  to  the  response  of  the  ionosphere  to  such  low- 
frequency  hydromagnetic  disturbances.  The  purpose  is  to  separate 
the  effects  of  the  medium  from  those  of  the  source  on  the  trans¬ 
mitted  disturbance. 
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SUMMARY 


In  this  report,  we  study  the  propagation  through  the  ionosphere 
of  sub -extremely -low  frequency  (<  0.1  cps)  hydromagnetic  waves.  We 
consider  the  case  of  plane  wave  propagation  in  the  vertical  direc¬ 
tion  in  the  presence  of  a  uniform  magnetic  field  which  is  also  in 
the  vertical  direction  (so-called  polar  propagation) .  The  purpose 
is  to  separate  out  the  effects  of  the  ionosphere  on  hycx omagnetic 
disturbances  propagating  through  it.  The  ionosphere  is  assumed  to 
extend  infinitely  far  in  the  upward  direction  and  to  terminate 
abruptly  80  km  above  a  perfectly  co’iducting  flat  earth.  The  charged 
particle  density  is  assumed  to  be  constant  throughout  the  ionosphere, 
and  the  neutral  particle  density  to  vary  exponentially  with  altitude. 
The  coupled  hydrodynamic  and  electromagnetic  equations  for  the 
medium  are  linearized,  and,  with  the  approximations  valid  for  low- 
frequency  waves,  the  resulting  system  is  shown  to  reduce  to  a 
second-order  wave  equation.  A  steady-state  solution  of  this  equa¬ 
tion  is  obtained,  from  which  transmission  and  reflection  coefficients 
cire  calculated.  It  is  found  that,  for  ionospheric  conditions 
corresponding  to  daytime  at  sunspot  maximum,  the  transmission  co¬ 
efficient  exhibits  a  pronounced  resonance  at  a  frequency  of  about 
0.013  cps.  Since  this  is  approximately  the  frequency  of  so-called 
Pc  oscillations,  it  is  suggested  that  the  resonance  just  described 
accounts  for  these  natural  micropulsations  in  the  geomagnetic  field. 
It  is  found  that  Vac.  resonance  shifts  to  higher  frequencies  at  night, 
also  in  agreement  with  the  data  on  natural  Pc  oscillations.  It  is 
shown  that  the  presence  of  an  almost  perfectly  conducting  ground 
(at  these  frequencies)  is  important  in  determining  the  location  of 
the  resonance.  Finally,  it  is  shown  that  disturbances  reaching  the 
ground,  in  the  polar  regions,  should  be  very  nearly  Icft-circularly 
polarized. 
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I.  INTRODUCTION 


A  subject  which  has  been  receiving  considerable  attention  for 
a  long  time  is  the  existence  of  sub -extremely -low  frequency  fluctua¬ 
tions  in  the  geomagnetic  field.  These  fluctuations,  which  have 
periods  of  the  order  of  seconds  to  tens  of  seconds,  are  associated 
with  many  naturally  occurring  geomagnetic  phenomena,  and  have  also 
been  generated  by  nuclear  explosions  at  high  altitude. 

The  experimental  data  on  natural  field  fluctuations  have  been 
summarized  by  Cambell  (1963)  .  Briefly,  the  data  show  regular 
oscillations  in  the  frequency  range  below  about  0.2  cps  appearing 
over  broad  sections  of  the  earth  with  related  phase  on  days  of  high 
solar-terrestrial  activity.  The  oscillations  appear  quite  sinu¬ 
soidal,  with  a  period  which  varies  slightly  with  season  and  con¬ 
siderably  with  the  time  of  day;  longer  period  oscillations  are 
observed  near  noon  whereas  shorter  periods  are  noted  at  night. 

Map^  (1959),  for  example,  found  two  daytime  bands  of  oscillations 
centered  at  70  seconds  and  20  seconds,  respectively,  and  a  nighttime 
band  with  an  average  peri^  if  8  seconds.  The  polarization  of  the 
magnetic  field  is  elliptical  in  the  horizontal  plane.  The  signal 
amplitude  shows  a  latitude  dependence,  dropping  off  toward  lower 
latitude  by  a  factor  of  about  \  per  20°  to  30°  lower  latitude. 

In  addition  to  the  data  on  natural  fluctuations,  much  data 
exists  on  the  fluctuations  generated  by  various  high-altitude 
nuclear  explosions.  This  data  shows  a  striking  similarity  between 
bomb-generated  and  natural  fluctuations  in  two  important  respects. 
In  the  first  place,  the  initial  phase  of  signals  generated  by  such 
explosions  arrives  essentially  simultaneously  at  stations  all  over 
the  world  (with  time  delays  after  detonation  which  depend  on  the 
altitude  of  the  explosion) .  Secondly,  the  signals  generated  by  a 
given  explosion'  are  essentially  the  same  in  shape  and  duration  at 
stations  all  over  the  world.  In  general,  one  may  say  that  the 
signals  resulting  from  low-yield  detonations,  such  as  the  Argus 
series,  are  quite  similar  in  most  respects  to  naturally  occurring 
signals. 
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One  further  point  of  Interest  is  the  power  spectra  associated 
with  oscillations  in  this  frequency  range.  Several  authors 
(Santirrocco  P,arkey .  1963;  Sn^^,  Pj^oyazek  and  Bp^jck^  1961; 
Can^bell^  1959;  and  tjorton  and  Hpffmann.  1962)  have  analyzed  the 
spectra  associated  with  natural  fluctuations,  while  tiodder  (1963) 
has  presented  a  spectral  analysis  of  signals  resulting  from  nuclear 
explosions  as  well  as  of  the  natural  background.  These  analyses 
show  one  or  more  large  peaks  in  the  frequency  range  of  interest. 

In  particular,  a  large  fraction  of  Hodder's  spectra,  both  bomb- 
induceci  and  naturally  occurring,  show  a  very  prominent  peak  at  about 
.02  cps. 

The  simultaneity  of  arrival  of  these  fluctuations  over  broad 
sections  of  the  earth  has  led  many  authors  to  conclude  that  the 
signals  received  at  aiious  stations  arise  from  a  common  ionospheric 
disturbance.  The  disturbance  is  generally  believed  to  propagate 
through  the  ionosphere  as  a  hydromagnetic  wave  which,  upon  emergence 
from  the  ionosphere,  propagates  as  an  electromagnetic  signal  to  the 
various  points  of  detection.  The  time  delay  between  detonation  and 
receipt  of  explosion-generated  signals  would  then  represent  the  trans¬ 
it  time  of  the  hydromagnetic  disturbance  through  the  ionosphere. 
Estimates  of  these  transit  times  are  in  reasonable  agreement  with  the 
observed  delays. 

The  point  which  concerns  us  in  this  paper  is  the  striking 
similarity  among  signals  in  this  frequency  range  resulting  from  a 
variety  of  natural  causes,  and  the  similarity  between  these  natural 
signals  and  those  which  are  produced  by  explosions.  These  similari¬ 
ties  have  suggested  to  many  authors  the  importance  of  the  ionosphere 
in  the  transmission  of  such  disturbances.  We  would  point  oat  here 
that  the  similarity  between  natural  fluctuations,  which  are  generated, 
for  the  most  part,  outside  the  ionosphere,  and  explosion-produced 
fluctuations,  which  are  generated  in  the  ionosphere  itself,  emphasizes 
the  particular  in^ortance  of  the  lower  ionosphere  (below  about  400  or 
500  km)  in  determining  the  characteristics  of  the  observed  signals. 

It  is  this  region  of  the  ionosphere  to  which  we  turn  our  attention. 
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The  low-frequency  transmission  characteristics  of  the  lower 
ionosphere  have  been  investigated  by  Francis  and  jtorplus  (1960)  and 
by  Field  (1963).  The  former,  making  certain  simplifying  assumptions, 
numerically  integrated  the  differential  equations  for  the  propaga¬ 
tion  of  hydromagnetic  disturbances  in  the  ionosphere,  using  empirical 
ionospheric  properties,  chosen  to  correspond  to  the  daytime  at 
sunspot  maximum.  Their  numerical  results  show  that  the  ionosphere 
is  essentially  transparent  for  angular  frequencies  somewhat  less 
than  1  radian/sec.  Fj^H  (1963)  used  empirical  ionospheric  proper¬ 
ties  to  calculate  an  altitude -dependent  complex  index  of  refraction 
for  the  ionosphere .  He  then  approximated  the  ionosphere  by  two 
layers,  the  upper  of  which  has  a  real  index  of  refraction  and  the 
lower  of  which  has  a  complex  index  of  refraction,  and  calculated 
the  resulting  transmission  coefficient.  Field’s  results  are  sub¬ 
stantially  in  agreement  with  those  of  Francis  -and  l^rplus  (1960)  . 
Neither  investigation,  however,  revealed  any  prominent  resonances 
in  the  transmission  coefficient. 

In  this  paper,  we  also  calculate  the  transmission  and  reflec¬ 
tion  coefficients  for  the  lower  ionosphere.  We  use  essentially  the 
same  model  as  Francis  Kaz^)li^  (1960),  but  with  some  additional 
simplifications  which  make  possible  an  analytic,  rather  than  a 
numerical,  solution.  The  advantages  of  an  analytic  solution  lie, 
of  course,  in  the  physical  insight  it  provides  and  in  the  ease  with 
which  the  effects  of  a  variation  in  ionospheric  parameters  can  be 
determined.  Even  with  the  necessary  simplifications,  the  important 
features  of  the  phenomenon  appear  in  such  a  model. 

We  start  by  incorporating  our  sin^lifying  assuiiq>tions  into 
the  basic  equations  of  motion  for  the  ionospheric  constituents. 

From  these  basic  equations  we  derive  a  wave  equation,  which  we  then 
examine  in  the  frequency  region  below  about  1  cps  for  the  case  of 
polar  propagation.  We  show  that,  for  the  polar  propagation  of  hydro- 
magnetic  waves  in  this  frequency  range,  the  ionosphere  consists 

I 

essentially  of  three  "regions":  a  lossless  "Alfven  region,"  above 

f 

about  350  km,  in  which  the  disturbance  propagates  as  a  pure  Alfven 
wave,  a  "Hall  region,"  below  about  130  km,  which  is  dominated  by 
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ion-neutral  particle  collisions  and  in  which  the  propagation 
characteristics  are  quite  different  from  the  Alfven  region,  and 
an  intervening  "transition  region"  which  smoothly  connects  the 
Alfven  and  Hall  regions.  (We  shall  sometimes  refer  to  the  Alfven 
region  as  the  "upper  ionosphere,"  and  to  the  transition  and  Hall 
regions  as  the  "lower  ionosphere.")  We  then  assume  that  the  Alfven 
speed  is  uniform  throughout  the  entire  (seni -inf inite)  ionosphere 
and  that  the  ion-neutral  particle  collision  frequency  (which  enters 
only  in  the  transition  region)  varies  exponentially  with  altitude. 

Wit.i  these  assumptions,  we  are  able,  without  the  further  introduction 
of  any  artificial  boundaries,  to  obtain  an  analytic  solution  of  the 
wave  equation.  The  solution  immediately  yields  the  desired  tranji- 
mission  and  reflection  coefficients.  The  transmission  coefficient 
is  found  to  exhibit  a  pronounced  resonance  at  a  frequency  of  .015 
cps  for  daytime  conditions,  with  an  upward  shift  in  frequency  for 
nighttime  conditions.  This  resonance  is  in  the  frequency  range  of 
so-called  Pc  oscillations,  which  are  therefore  apparently  accounted 
for  by  the  filtering  action  of  the  lower  ionosphere.  (This  resonance 
does  not  appear  in  the  work  of  Francis  ^rplus,  (1960)  ,  whose 
calculations  extend  only  to  a  frequency  of  1  rad/sec  at  the  low  and, 
nor  in  the  work  of  Field  (1963)  who  was  concerned  only  with  an 
average  transmission  coefficient  rather  than  in  its  detailed  dependence 
on  frequency.) 

9 

It  should  be  mentioned  that  our  assumption  of  a  constant  Alfven 
speed  throughout  a  serai-infinite  ionosphere  treats  the  upper 
ionosphere  as  a  lossless,  non -dispersive  medium  (or  transmission 
line).  It  is  known,  however,  that  the  Alfven  speed  actually  in¬ 
creases  with  altitude  in  the  region  between  about  200  km  and  2000  km. 
This  should  be  of  little  consequence  for  the  initial  phase  of  bomb- 
induced  signals  which  are  generated  at  altitudes  of  a  few  hundred 
kilometers,  but  may  be  expected  to  modify  signals  which  are  generated 
at  altitudes  of  2000  km  or  more.  The  transmission  properties  of  a 
medium  with  varying  Alfven  speed  have  been  investigated  by  several 
authors,  and  Wat  ana  be  (1962)  have  considered  the  case  of  a 

collisionless  upper  ionosphere  in  which  the  Alfven  speed  increases 
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exponentially  with  altitude  to  a  maximum  at  2000  km.  They  approxi¬ 
mated  the  medium  by  six  uniform  Alfven  layers  of  equal  thickness, 
the  lowest  of  which  is  grounded,  and  numerically  calculated  a  trans¬ 
mission  coefficient,  ^iel^d  (1964-)  has  obtained  an  analytic  solution 
to  the  same  problem.  These  authors  find  that  a  filtering  action 
takes  place  between  the  region  of  maximum  velocity  at  around  2000  km 
and  th-'-  bottom  of  the  upper  ionosphere.  Various  resonances  appear, 
the  lowest  of  which  occurs  at  a  frequency  of  about  .1  cps,  or  about 
an  order  of  magnitude  higher  than  the  frequency  of  the  resonance 
associated  with  the  lower  ionosphere.  Prince  and  Bostick  (1964) 
have  carried  out  numerical  calculations  using  a  ir.nre  realistic 
ionosphere,  with  various  collisional  effects  taken  into  account. 

They  find  similar  resonances  in  the  transmission  coefficient.  These 
resonances  in  the  transmission  coefficient  of  the  upper  ionosphere 
are  thought  to  be  responidble  for  ’’pearl"  type  oscillations. 

For  disturbances  originating  at  heights  of  2000  km  or  more  (and 
in  the  frequency  range  below  about  1  cps)  ,  the  ionosphere  below  2000 
km  therefore  appears  very  much  like  two  successive  filters.  The 
upper  filter  has  a  transmission  resonance  at  about  .1  cps,  and  the 
lower  filter  a  transmission  resonance  at  about  .015  cps.  The  com¬ 
bined  transmission  coefficient  is  essentially  the  product  of  the 
separate  transmission  coefficients. 
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II.  WAVE  EQUATION  IN  THE  IONOSPHERE 


We  start  by  making  the  same  simplifying  assumptions  as  Francis 
and  ^yplus_  (1960),  viz., 

(1)  the  static  magnetic  field  is  uniform, 

(2)  the  pressure  is  so  low  that  hydrostatic  restoring  forces 
and  gravitational  forces  are  small  compared  to  electro¬ 
static  restoring  forces, 

(3)  the  fractional  ionization  is  so  small  that  collisions 
between  electrons  and  ions  may  be  neglected, 

(4)  the  fractional  ionization  is  so  small  that  collisions  of 
neutrals  with  charged  particles  may  be  neglected  (i.e., 
the  neutrals  do  not  participate  in  the  mass  motion) ,  and 

(5)  the  ionosphere  is  still  and  flat,  and  the  wave  normal  is 
along  the  ionospheric  density  gradient. 

In  addition  to  these  assun^tions,  we  further  assume  that 

(6)  the  static  field  is  uniform  and  in  the  vertical  direction, 
which  restricts  the  applicability  of  our  model,  la  detail  at 
least,  to  the  polar  and  auroral  regions.  These  basic  assuiiq>tions 
will  now  be  incorporated  into  the  equations  of  motion  for  the 
ionospheric  constituents. 

We  start  with  the  following  transport  equations  for  the  electrons 
and  (singly -charged)  ions: 


V  X  B 


at  ^i 


V  v; 


-  —  (E  + 
”i 


(la) 


(lb) 


where  and  are  the  collision  frequencies  of  electrons  with 
neutrals  and  ions  with  neutrals,  respectively,  and  the  other  symbols 
have  their  usual  meanings.  (We  shall  use  cgs  units  throughout.)  If 
we  take  the  positive  z-direction  to  be  vertically  downward  and 
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restrlct  our  attention  to  plane  transverse  waves  of  small  amplitude. 


we  may  write 


V  -  (Vj^,  v^,  0) 


B  -  (b^,  b^,  B^) 


E  -  (E^,  E^,  0)  (2c) 

where  the  v,  b,  and  E  are  assumed  to  be  small  quantities.  Linearizing 
the  transfer  equations,  we  obtain 


“  +  u)  V  X  e_  +  j  V 

3t  e  e  3  ‘'e  e 


—  E 
m 

e 


-rr  -  uo,  V.  X  e.  +•  V,  V.  -  —  E 
at  i  1  3  i  1 


where  e,  ie  a  unit  vector  in  the  positive  z -direction,  and  o)  ,  »  — 

3  e,l  jC 

is  the  cyclotron  frequency  of  the  particular  ionospheric  component. 

With  the  assumption  of  a  static  field  in  the  vertical  direction,  an 

essential  simplification  of  Che  equations  results  if  we  introduce 

the  complex  variables 


V  -  Vi  -  1  Vz 


E  -  E^  -  i  E^ 


In  terms  of  these  variables ,  Eqs .  (3a)  and  (3b)  take  the  form 


(fig  ^  i  Vg 


(^1  -  i  (Bj) 


—  E 
®l 


(5b) 
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where  fj  is  the  operator 

n  H  +  u  (6) 

for  the  component  in  question. 

Remembering  that  the  current  density  is  defined  as 

J  “  n  e(v^  -  v^)  (7) 

where  n  is  the  number  density  of  electrons  and  ions  (assumed  equal), 
we  may  combine  Eqs.  (5a)  and  (Sb)  to  obtain 

(fi^  +  iu)g)(Q.  -  i  =  «  [  ~  +  iuJg)  +  ~  E  (8) 

L  i  “g  J 

where 

J  -  -  I  jj.  (9) 

We  note,  in  passing,  that  Eq.  (8)  provides  a  very  simple 
derivation  of  two  of  the  diagonal  components,  and  of  the 
conductivity  tensor.  The  equation  defining  the  conductivity  tensor 
reduces,  in  the  present  case,  to 

j  =  E  -»■  02  E  X  Cj  (lO) 

which  may  also  be  written 


j  =  (a^  +  i  aj)  E.  (11) 

*  i  'jut 

Assuming  a  time  dependence  of  the  form  e  ^.n  Eq .  (8)  ,  and 
conq>aring  Eqs.  (8)  and  (11),  we  obtain 


n  e 


( 


-  i» 


^i 


-  iu; 


".[(Oi- 


i»)' 


2-1 

*3 


*Bi[(nj^”i«>)  + 


.2, 


(12a) 


0). 
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2  r 

n  e  •<  - 


U) 


x^1 


2  2, 

tn.  rCu.  -i x)  x.  1 

1-  i  i"* 


(12b) 


in  agreement  with  the  usual  formulas  for  these  conductivities. 

Finally,  to  derive  a  wave  equation  for  E,  we  combine  Eq .  (8) 
with  Eq.  (7)  and 


E  = 


=  in  21 

2  ^t 
c 


(13) 


to  obtain 


+  iU)^)(fi. 


ia-p 


E 


4tt  ne 


m  m.c 
e  1 


(in  Q 
e  e 


la.O.) 
1  1 


2E 

at 


(14) 


where  we  have  used  m  uu  =  m.x. .  The  general  wave  equation  is  evi- 
dently  a  complicated  fourth-order  equation  and  we  shall  not  attempt 
to  discuss  its  general  properties  at  this  time.  We  do  note,  how¬ 
ever,  that  in  the  limit  of  very  strong  magnetic  field  and  no 
collisions,  i.e.,  uu  ,  x.  -•  »  and  u  ,  v.  -  0,  Eq .  (14)  reduces  to 

I 

the  usual  equation  for  Alfven  waves 

V^E=\^.  (15) 

V  at 
a 


where  v  = 


‘♦Tt 


B 

o  * 

— 7 — - ^  is  the  square  of  the  Alfven  speed. 

n(m.+  m  ) 

1  e 


Fortunately,  for  the  frequency  regime  under  consideration, 

Eq .  (14)  simplifies  considerably.  Much  of  the  simplification  re¬ 
sults  from  the  properties  of  the  ionosphere  in  the  region  of 
interest  (80  km  -  500  km),  viz., 

(a)  X  ~  10^  rad/sec, 

'  c  * 
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(b)  10^  rad/sec, 

2 

(c)  0;^  ~  10  rad/sec, 

(d)  10  ^  <  10^  rad/sec  - 

(Sae  Francis  and  Karp^s  (1960)  for  a  tabulation  of  ionospheric  proper¬ 
ties).  In  view  of  these  relationships,  and  the  fact  that  we  are 
interested  in  frequencies  a;  «  (i.e.,  «  •JU^)  »  the  following 

approximations  become  valid: 


fl  +  ix  ix 

e  e  c 

(16a) 

-  ‘“i  *  »i  -  • 

”e  ^ 

Furthermore,  since  >  ( — )  o  , 

1  m.  e 

i 

(16b) 

m  n  m.  n.  ft-  m,  P.. 

e  e  11  i  i 

(16c) 

L 

to  terms  of  relative  order  (m  /m.)^. 

With  these 

approximations , 

Eq.  (14)  becomes 

2  1 
(u.  -  ix.)  V  E  *  -  ix.  “T 
1  1  1  2 

n  ^ 

“i  St  * 

(17) 

V 

a 

a  somewhat  more  tractable  second-order  equation. 
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Before  attempting  a  solution  of  Eq.  (17),  it  is  worth  noting 
that  for  hydromagnetic  waves  in  the  fre.juency  regime  below  about 
1  cps ,  the  ionosphere  consists  of  three  essentially  distinct  regions, 
which  are  characterized  by  the  relative  magnitudes  of  the  three 
frequencies  and  x  (the  angular  frequency  of  the  disturbance) . 

In  order  of  increasing  altitude,  these  regions  are: 

(a)  Hall  Region  -  This  is  defined  to  be  the  region  in  which 
the  Inequality  x  «  x\  <  u.  is  satisfied,  i.e.,  where  the  collision 
frequency  is  the  dominant  frequency.  It  lies,  therefore,  between 
about  80  km  and  130  km.  The  boundaries  of  this  region  are  inde¬ 
pendent  of  x,  provided  that  x'  «  Here,  the  Hall  effect  pre¬ 

dominates;  in  terms  of  the  conductivities  defined  by  Eqs.  (11a)  and 
(11b),  02  »  In  this  region,  Eq.  (17)  reduces  to 


E 


•  1  dE 

^  2  At  * 

V 

a 


(18) 


i.e.,  to  a  Schrodinger-like  equation.  It  is  important  to  note  that, 
though  this  region  is  collision-d.  minated,  the  wave  equation  never¬ 
theless  does  not  contain  the  collision  frequency  explicitly.  The 
implicit  inq)ortance  of  collisions  is,  however,  reflected  in  the  form 
of  the  equation. 

(b)  Transition  Region  -  This  region  is  defined  by  the  inequality 

u)  <  T3.  <  uu. .  It  lies  between  about  130  km  and  353  km.  The  location 
1  1 

of  the  upper  boundary  of  this  region  depends  on  the  frequency  uu; 
the  value  given  is  appropriate  to  a  frequency  x  st  .3  rad.^sec.  Here 
no  further  simplification  of  the  wave  equation  is  possible;  collisions 
thus  play  an  explicit  role  in  this  region. 

(c)  Alfven  Region  -  This  is  the  region  in  which  the  inequality 
<  X  «  'jUj^  is  satisfied.  It  is  therefore  the  region  above  about 

330  km,  and  in  our  model  is  assumed  to  be  semi-infinite  in  extent. 

In  this  region,  Eq.  (17)  reduces  to  Eq.  (14),  the  usual  equation 
for  Alfven  waves.  The  wave  equation,  in  this  region,  is  independent 
of  the  collision  frequency. 


actual  functional  form  of  the  collision  frequency  enters  only  in 
the  transition  region.  To  obtain  an  analytic  solution  to  the  wave 
equation,  we  must,  of  course,  approximate  the  ionospheric  collision 
frequency  by  an  analytic  function.  To  obtain  an  analytic  approxi¬ 
mation  which  is  valid  over  a  large  region  of  the  ionosphere,  and 
which  at  the  same  time  permits  an  analytic  solution  of  the  wave 
equation,  does  not  seem  possible..  We  see,  however,  that  this  is 
fortunately  not  necessary;  we  are  presented  with  the  far  simpler 
task  of  constructing  an  expression  which  is  a  good  approximation 
only  over  a  relatively  limited  portion  of  the  ionosphere. 

Finally,  let  us  emphasize  that  we  have  been  speaking  of  "regions" 
in  a  qualitative  rather  than  in  a  mathematical  sense.  In  the  model 
which  we  employ  below,  it  will  not  be  assumed  that  these  regions 
are  separated  by  artificial  sharp  boundaries,  with  a  separate  wave 
equation  appropriate  to  each.  The  single  wave  equation,  Eq.  (17), 
will  govern  the  propagation  of  the  disturbance  throughout  the  entire 
ionosphere.  The  essence  of  the  above  discussion  is  that  the 
character  of  the  wave  equation  changes  (continuously,  of  course) 
from  one  region  to  the  next,  and  this  property  will  obviously  be 
reflected  in  the  solution  to  that  equation. 


-13- 


SGwaIGK  Or  THE  WAVE  EQuivii-uw  tuts,  a  lu/uau  lUNUstnejm, 


In  this  section  we  shall  construct  the  steady-state  solution  to 
Eq.  (17)  for  a  model  ionosphere  which  will  be  described  below.  If 
we  introduce  a  time  dependence  of  the  form  e  Eq.  (17)  becomes 

2  2 

~  n^  (u);z)  E  *  0  (19) 

dz  c 


where 


n^(uu;z) 


^2  (1  -  iuj/uj^) 

^2  u)  (1  -  icoj^/oj^)* 

a 


(20) 


In  Eq.  (19),  E  now  denotes  E(z,(u),  the  Fourier  transform  of  E(z,t). 
With  a  time  dependence  of  the  form  e  and  the  definition  of 
E(z,t)  by  Eq.  (4b),  positive  values  of  o)  correspond  to  left-circular ly 
polarized  waves,  which  we  shall  call  the  L-modc,  and  negative  values 
of  (u  correspond  to  right-circular ly  polarized  waves,  which  we  shall 
call  the  R-mode.  (In  the  case  of  polar  propagation  which  we  are 
considering,  these  two  modes  are  decoupled  from  each  other.) 

The  more  important  properties  of  our  three  regions  can  already 
be  deduced  from  Eqs.(19)  and  (20).  In  the  Alfven  region,  U)/i>£»  1 
and  u-j/vj  »  1*  so  that  the  index  of  refraction  n(u);z)  becomes 

n(u);z)  -  "n  *  ^  •  (21) 

a 


This  is  the  well-known  result  that,  for  both  modes  of  propagation, 

I 

the  Alfven  region  is  a  lossless  and  non-dispersive  medium. 

In  the  Hall  region,  uu/uj^  «  1  and  «  1»  the  index 

of  refraction  becomes 


«i  k 

n(u);z)  njj(u))  -  (~)  . 


(22) 
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This  !aedi'«2  thetjfcre  appears  quite  ulffcreni.  to  the  tvfo  modes  of 
propagation.  For  the  L-mode  (ud  >  0) ,  the  index  of  refraction  is  real, 
and  this  mode  therefore  propagates  through  the  Hall  region  without 
attenuation,  but  obviously  with  dispersion.  For  the  R-mode  (ou  <  0) , 
the  situation  is  quite  different;  the  index  is  pure  imaginary,  and 
this  mode  is  consequently  evanescent.  Because  of  this,  the  trans¬ 
mission  coefficient  for  the  R-mode  drops  off  very  rapidly  with  fre¬ 
quency,  and  it  is  mainly  the  L-mode  which  gets  through  the  ionosphere. 
The  disturbances  reaching  the  ground  should  thus  be  circularly 
polarized.  (We  must  remember,  of  course,  that  we  are  treating  the 
case  of  polar  propagation;  in  the  higher  latitude,  but  non-polar, 
regions^  we  would  expect  the  same  mechanism  to  produce  eliptical 
polarization,  in  agreement  with  observation.) 

In  the  transition  region,  no  similar  simplification  of  the 
index  of  refraction  takes  place.  The  index  for  both  modes  is  complex, 
so  that  in  this  region  each  mode  undergoes  both  attenuation  and 
dispersion. 

2 

The  quantities  and  v  in  Eq.  (19)  are,  of  course,  all 

X  X  d 

functions  of  the  height  z.  The  Alfven  speed  v  depends  on  the 
density  of  ions,  which  varies  with  altitude  (and  time  of  day).  How¬ 
ever,  as  we  have  pointed  out,  we  are  concerned  only  with  the  trans¬ 
mission  properties  of  the  lower  ionosphere  (i.e.,  below  about  350  km). 
Over  this  region  the  ion  density  is  relatively  constant  (except  at 
the  very  bottom,  where  we  will  replace  the  very  rapid  actual  varia¬ 
tion  with  a  discontinuity) .  We  will  assume  that  the  Alfven  speed  is 
constant  throughout  the  whole  ionosphere,  with  a  value  appropriate 
to  the  lower  ionosphere.  As  we  have  stated  previously,  this  approxi¬ 
mation  should  have  little  effect  on  the  calculated  properties  of  the 
lower  ionosphere.  The  same  remarks  apply  to  the  ion  cyclotron  fre¬ 
quency  which  is  also  assumed  to  be  constant. 

It  now  remains  to  approximate  the  ion -neutral  particle  collision 
frequency,  For  reasons  discussed  above,  it  is  sufficient  to  have 

an  analytic  approximation  which  is  accurate  only  over  the  transition 
region.  Over  this  region,  the  collision  frequency  is  approximated 
rather  well  by  an  exponential;  we  therefore  take  to  have  the  form 
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i 


(V. 


(23) 


We  have  taken  the  plane  z  >■  0  to  be  at  that  altitude  at  which 

I 

or  In  other  words,  at  the  "boundary"  between  the  Alfven  and  transi¬ 
tion  regions,  which  occurs  at  an  altitude  of  about  130  km.  (The 
positive  z-direction  has  been  taken  vertically  downward,  as  before.) 
For  the  time  being,  the  collision  frequency  scale-height  Q  will  be 
left  as  a  parameter,  with  an  appropriate  numerical  value  to  be 
assigned  later.  This  will  enable  us  to  examine  the  role  played  by 
the  transition  region  in  determining  the  transmission  coefficient. 

It  is  convenient  at  this  point  to  work  with  dimensionless  quanti¬ 
ties,  which  we  shall  denote  by  a  tilde.  To  do  so,  we  introduce  the 
following  units : 

unit  of  frequency  «  u)^ 

V 


unit  of  length 


In  these  units,  the  frequency  range  of  interest  is  u)  «  1.  Intro¬ 
ducing  these  units  and  combining  Eqs.  (19),  (20),  and  (23),  we  obtain 


^  *  t(S;  i/C)  e  -  0 

dz 


where 


f(w;  r/c) 


u)  -  itt)  e  ^ 

* 

1  -  ie-*/C 


(24) 


(25) 


The  propagation  problem  actually  extends  down  to  the  surface  of 
the  earth,  which  we  treat  as  a  perfect  conductor  at  the  low  frequencies 
of  interest  to  us.  At  the  base  of  the  ionosphere,  the  variation  in 
composition  is  very  rapid  compared  with  the  wavelength  in  its  vicinity; 
we  therefore  replace  the  actual  variation  in  composition  by  a  dis¬ 
continuity  at  z  •  z^.  With  our  choice  of  coordinate  axes,  the 
atmosphere  is  therefore  divided  in  the  following  way  into  two 
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atmospheric  regions  (the  upper  of  which  we  have  already  subdivided 
into  ionospheric  regions) ; 

2  <  Zj^;  ionosphere 
(2*0);  altitude  at  which 

2*2^;  discontinuous  transition  to  neutral  atmosphere 
Zj^  <  2  <  Zj^  +  h;  neutral  atmosphere 
2  *  Zj^  +  h;  surface  of  the  earth. 

The  quantity  h  is  the  height  of  the  neutral  atmosphere. 

Since  the  earth  is  a.ssumed  to  be  a  perfect  conductor,  the  field 
at  the  surface  of  the  earth  must  vanish.  In  the  neutral  atmosphere, 
the  field  must  therefore  have  the  foim 

E  *  ^  (26) 

(S'.  <  r  <  2  ^  +  to 


where  A^(a5)  is  some  unknown  function  of  oj.  The  boundary  condition 
we  iiiq>ose  at  the  base  of  the  ionosphere  is  that  the  field  and  its 
derivative  are  continuous,  or  that 


1 


-ir  cotir  tr 

o  o 


(27) 


where  the  approximation  of  the  cot  follows  from  the  fact  that  the 
wavelength  in  air  is  extremely  long. 

We  next  turn  our  attention  to  the  solution  of  the  wave  equation 
in  the  ionosphere.  From  Eq.  (25)  it  is  apparent  that  for  z*/^  »  1 
(i.e.,  in  the  Hall  region),  f((B;  iu.  There  must  consequently 

exist  two  linearly  independent  solutions,  and  W^,  cf  Eq.  (24) 
aach  that 


W 


1 


W 


2  -•  uc  . 


(28) 
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We  take  our  general  solution  of  Eq.  (24)  to  have  the  form 


E  *  A^(^S>(W^  +  oWp 


(29) 


where  A^^Cuu)  and  a  are  two,  as  yet  unknown,  functions  of  w.  Since  the 
asymptotic  form  given  by  Eq .  (28)  is  a  good  approximation  in  the  Hall 
region,  we  have 


E 


A^(yE) 


(e 


+  eye 


-lU) 


(30) 


Using  Eq .  (30)  to  evaluate  the  left-hand  side  of  Eq.  (27)  and  solving 
for  a,  we  obtain 


Of 


Cl-h  i^  ^  Tv) 

(1  -  iuj  ^  hi) 


(31) 


which,  substituted  back  into  Eq.  (30)  gives,  for  the  electric  field 
at  the  base  of  the  ionosphere. 


E(zj^;w) 


1  ^ 
1:1  .l.nj-  e  1 
(1  -  iTj;  ^  10 


(32) 


To  obtain  a  transmission  and  reflection  coefficient,  ve  must 
examine  our  solution,  Eq.  (30),  in  the  Alfven  region  (z/^  «  -  1). 

In  the  Alfven  region,  Eq.  (25)  becomes  f(u);  ps  u)  ;  consequently, 

any  solution  of  Eq.  (24)  teconies  some  linear  combination,  in  the 

I  jp 

Alfven  region,  of  the  two  linearly  independent  solutions  e  and 
e  Ve  may  therefore  write 


lU)  z  .  , 
e  4  b, 


lU)  z  .  , 
e  4  b, 


-iu3  z 
e 

-lU)  z 
e 


(Z*  -  -  oo)  . 


(33) 
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Ccmbining  Eqs.  (33)  and  (29),  we  obtain  for  the  field  in  the  Alfven 
region 


Efi-AjCu^  (a^  +0B2)(e  ■♦■re  )  » 

(?  -  -  oo) 


(34) 


where 


r 


a. 


+  cyb 
+  jya 


2 

2 


(35) 


We  shall  refer  to  |r|  as  the  anq>litude  reflection  coefficient;  the 
square  of  this  coefficient  is  the  usual  energy  reflection  coefficient. 

ProK  Eqs.  (30)  and  (34),  we  can  calculate  the  ratio  of  the 
electric  field  at  the  base  of  the  ionosphere  to  the  incoming -wave 
part  of  the  electric  field  in  the  Alfven  region.  This  quantity, 
which  we  call  the  aii^illtude  transmission  coefficient  for  the  electric 
field,  and  which  we  denote  by  t_,  is  given  by 


E(z^,^ 

Aj(uj)  (a^^harap 


.  ltt>  z 
2i^  ^  TT  e  ^ 

(1-itB  ^  K)(aj^-t*  Qre2) 


(36) 


In  practice,  we  are  interested  in  the  magnetic  field  at  the 
surface  of  the  earth  which,  at  these  extremely  large  wavelengths  in 
air,  is  essentially  the  saate  as  the  magnetic  field  at  Che  base  of 
the  ionosphere.  Combining  the  Induction  equation  with  Eq.  (27),  we 
obtain 


(37) 


where  B  •  B  -  i  B  .  On  the  other  hand,  the  induction  equation 
a  y 

applied  to  Eq.  (34)  gives 
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B  »5  icA. (u))(a 


+002) 


lu;  z 


-re 


-10)  z. 


(z - ®) 


fto  that,  In  analogy  with  Eq.  (36),  the  amplitude  transmission  co¬ 
efficient  for  the  magnetic  field,  which  we  denote  by  t^,  is  given 
by 


iou  z. 


2  e 


^B 


ciA^(^(aj^4-  ara2) 


0)  ^(1-iu)  ^  K)  (aj-‘-Qra2) 


(38) 


It  remains  now  only  to  determine  the  quantities  a^,  bj^, 
and  b2.  These  are  obtained  by  determining  the  solutions,  and 
W2,  of  Eq.  (24)  which  have  the  form  prescribed  by  Eq .  (28)  and 
examining  them  in  the  limit  z'  where  they  must  assume  the 

form  given  by  Eq.  (33).  In  order  to  avoid  introducing  any  un¬ 
necessary  mathematical  detail  at  this  point,  we  shall  defer  this 
procedure  to  the  Appendix  and  merely  present  the  results  in  the 
next  section. 
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The  numerical  results  presented  in  this  section  will  be  based 
on  the  following  values  of  the  relevant  parameters,  chosen  to  coin¬ 
cide  with  the  ionosphere  of  Francis  K.^rp lus  (1960)  : 

h  “  80  km 
*  50  km 

(1)^  (unit  of  f'-e^uency)  »  175  rad/sec 

^a  1.2 

—  (unit  of  length)  =  — r 

•"i  N.^ 

1 

”6 

“  ion  density  in  (no./cc)  x  10 

We  recall  that  the  cyclotron  frequency  u)^  is  considered  constant 
throughout  the  ionosphere ;  the  numerical  value  chosen  for  is 
that  frequency  at  which  in  the  actual  ionosphere,  which,  by 

Eq.  (23),  also  defines  the  plane  z  =  0.  The  numeric  in  the  unit  of 
length  is  based  on  an  ion  atomic  weight  of  26,  which  is  also  appro¬ 
priate  to  the  region  around  z  »  0,  i.e.,  around  130  km. 

To  obtain  the  transmission  and  reflection  coefficients,  we  must 
introduce  the  expressions  obtained  in  the  Appendix  for  aj^,  32*  bj^, 
and  b-  into  Eqs .  (38)  and  (35).  The  results  are  most  conveniently 

^  V 

expressed  in  terms  of  a  function  g((a  ,u0  defined  by 


gCm 


UJ  ^)  r.  #TI  h 

— ^ ' — 2  (1  -  lUJ  ^  exp(2  t  ■ 


rr(i-iT^  ^ 


■J'  U)  -  iz"^^ 


where,  as  usual,  F  denotes  the  Gamma -function.  Combining  the  results 
of  the  Appendix  with  Eqs.  (35)  and  (38),  respectively,  we  obtain 

r  »  -  II  -  -^)  gC(P  ^  Cl  ^  h 

(1  4*  U)  g(uj  4-  (1  -  UB  g(-^  ^,0)) 
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ra  -  2i?^5[a  *T,h  gGs  *',3)  +  a  -  s  h  8(-s  ^35]  ’ 


(41) 


(As  far  as  the  functional  form  of  g  is  concerned,  cu  and  u)  are 
created  as  independent.) 

Before  obtaining  numerical  results  for  r  and  t^,  it  will  be 
instructive  to  examine  the  limiting  case  *  0.  As  can  be  seen 
from  Eq.  '23),  this  case  corresponds  to  a  collision  frequency 


1)^-0  z  <  0 

(42) 

“  z  >  0  , 

and  therefore  to  a  medium  consisting  of  two  uniform  layers,  an 
Alfven  region  and  a  Hall  region,  (in  which,  despite  the  infinite 
collision  frequency,  it  is  still  assumed  that  the  neutral  particles 
do  not  participate  in  the  mass  motion),  separated  by  a  sharp  boundary 
at  z  »  0.  (This  limiting  case  can,  of  course,  be  obtained  much 
more  simply  by  direct  calculation.)  With  ^  *  0,  all  of  the  Gamma- 
functions  in  the  above  expressions  are  unity,  and  considerable 
simplification  results.  It  is  not  difficult  to  show,  for  example, 
that 


|r|  =  1  (43) 

for  both  the  L-mode  and  the  R-mode,  although  the  physical  signifi¬ 
cance  of  this  result  is  quite  different  for  the  two  modes.  For 
the  L-mode,  this  is  an  expression  of  the  fact  that  both  regions  are 
lossless  (real  index  of  refraction  in  both  regions) ,  and  the 
perfectly  conducting  earth  reflects  all  of  tbe  incident  energy.  On 
the  other  hand,  the  R-mode,  as  we  have  pointed  out  previously,  has 
a  pure  imaginary  index  of  refraction  in  the  Hall  region,  and  is  con¬ 
sequently  evanescent  in  this  region.  This,  coupled  with  the  boundary 
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conditlon  which  we  have  imposed  at  the  surface  of  the  earth,  makes 
this  evanescent  region  look  infinite  to  the  incoming  wave  (although 
with  an  index  of  refraction  somewhat  different  from  that  of  the 
actual  Hall  region),  and  consequently  leads  to  total  reflection. 

The  transmission  coefficient,  in  this  limiting  case,  is  given 
by 


'B 


-iz.u) 


iz.uj 


[(1+^^)(1  -  itt)^!0  e  ^  +  (1  -  ^  S(1  h  ^  TO  e  ^  ] 


(44) 


This  coefficient  will  exhibit  resonances  at  frequencies  for  which 
the  denominator  of  Eq.  (44)  becomes  small.  It  is  a  well-known 
result  that  the  nature  of  such  resonances  is  intimately  related  to 
the  location  in  the  conqjlex  u^plane  of  the  poles  of  the  trans¬ 
mission  coefficient.  If  a  pole  lies  close  to  the  real  m-axis 
(and  we  shall  show  shortly  that  this  is  indeed  the  case  for  the 
poles  of  Eq.  (44)  in  the  frequency  range  of  interest)  it  represents 
a  resonance,  the  real  part  of  the  pole  giving  the  position  of  the 
resonance  and  the  ioiaginary  part  giving  its  width. 

Now,  the  poles  of  Eq.  (44)  are  the  roots  of  the  equation 


(1  +  w  ^)(1  -  m  *  TO  e  +  (1  -  m  ^)(1  +  im  ^  TO  e^ 


0. 


(45) 


We  are  interested  in  the  frequency  range  ud  ^  «  1,  where  Eq,  (45) 
becomes,  approximately. 


cot  (m  ^  z'p  -  m  ^  T»  *  0. 


(46) 


It  is  not  difficult  to  show  that  the  roots,  u)^,  of  this  approximate 
equation  are  pure  real  and,  consequently,  that  the  low-frequency 
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poles  of  £q.  (44)  lie  very  close  to  the  real  (v-axis,  as  was  scared 
above.  These  poles  thus  represent  very  sharp  resonances,  as  can  be 
seen  in  Pig.  1  where  |tg|,  for  the  L-mode,  is  plotted  against  uj  for 
the  case  C  *  ®  (with  z.  »  50  km  and  h  ■  80  kro) ,  for  ion  densities 

5  5  3 

of  5  X  10  /cc,  10  /cc,  and  5  x  10  /cc.  The  three  values  chosen  for 
the  ion  density  are  representative,  in  order  of  increasing  density, 
of  conditions  in  the  daytime  at  sunspot  maximum,  in  the  daytioie  at 
sunspot  minimum,  and  in  the  nighttime,  respectively. 

A  simple  physical  interpretation  of  these  resonances  exists  in 
the  limiting  case  where  we  ground  the  ionosphere  at  its  base,  i.e., 
for  the  case  h  *  0.  In  this  case,  the  roots  of  Eq.  (46)  are  given  by 


n 


(2n-l)  J  , 

n  *  1,  2,  ... 


(47) 


where  is  the  wavelength  in  the  Hall  region.  The  lowest  root 

therefore  corresponds  to  the  frequency  at  which  the  Hall  region  is 

a  quarter-wave  plate.  However,  it  should  be  noted  that,  since  the 

term  ud  ^  ti  in  Eq.  (46)  is  of  order  unity  even  for  the  lot/est  of  these 

roots,  the  inclusion  of  a  neutral  atmosphere  between  the  base  of  the 

ionosphere  and  ground  plays  a  very  iirportant  role  in  determining  the 

positions  of  the  resonant  frequencies.  For  example,  for  the  case 

of  an  ion  density  of  5  x  lO^/cc,  the  lowest  frequency  determined  by 

-3 

Eq.  (47)  is  a  *  3  X  10  ,  whereas  '.be  lowest  root  of  Eq.  (46)  is 

-4 

actually  at  w  »  6  x  10  ,  a  factor  of  5  lower.  Thus,  the  resonances 

are  shifted  toward  significantly  lowex  frequencies  by  the  inclusion 
of  a  region  of  neutral  atmosphere.  This  emphasizes  the  importance 
of  including  such  a  region  in  any  calculation  of  ionospheric  trans¬ 
mission  coefficients. 

Before  leaving  this  two-layer  limit,  we  note  from  Fig.  1  that 
the  transmission  resonances  show  an  upward  shift  in  frequency  from 
day  to  night,  and  from  sunspot  maximum  to  sunspot  minimum,  as  is  ob¬ 
served  for  Pc  oscillations.  This  result  persists  even  when  a  transi¬ 
tion  region  is  introduced.  This  is  best  understood  by  rewriting  Eq. 
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(46)  in  terms  of  the  dimensional  lengths  z,  and  h.  recalline  that 
the  unit  of  length  is  v  /u>  «  c  N  ,  where  c  is  a  constant  which 

A  1  X  1  1 

depends  only  on  the  ion  mass.  In  terms  of  these  quantities,  Eq . 
(46)  becomes 


cot  (u)  ^  N  ^  c.z.)  - 
'  ill 


U) 


N,  ‘ 
1 


0. 


(48) 


It  is  clear  that  the  roots  of  this  equation  can  be  written 


u)  ^  N>  =  d 
n  1  n 


(49) 


where  the  d^  are  constants  independent  of  ionospheric  conditions. 

We  thus  find  that  each  resonant  frequency  varies  inversely  with  ion 
density.  This  result  holds  strictly  in  the  two-layer  limit  only  for 
uj^<  since  Eq.  (46)  was  derived  from  Eq.  (45)  on  the  basis  of 
that  assun^tion.  With  the  same  limitation  on  u),  this  same  result 
can  be  shown  to  apply  in  the  general  case  to  the  roots  of  Eq.  (41), 
as  well. 

The  two-layer  limit  considered  above  has  two  major  inadequacies 
it  introduces  an  artificial  sharp  boundary  between  the  Alfven  and 
Hall  regions,  which  may  give  rise  to  spurious  resonances,  and  it 
fails  to  include  the  effects  of  collisional  losses,  the  very  region 
in  which  such  losses  occur  having  been  eliminated  in  this  limit. 

To  include  a  transition  region,  we  must  calculate  the  transmission 
and  reflection  coefficients  for  an  appropriate  value  of  the  scale 
height,  Q.  There  seems  to  be  considerable  disagreement,  however, 
as  to  what  the  appropriate  ionospheric  scale  height  is .  From  the 
ionosphere  of  Francis  and  Radius  (1960),  for  example,  a  scale 
height  of  14  km  would  appear  appropriate,  whereas  the  ionosphere  of 
l^ince  Md  Bostick  (1964)  puts  the  scale  hei^t  closer  to  10  km. 
Because  of  this  uncertainty  in  the  best  value  of  we  shall  present 

the  results  for  both  ^  =  10  km  and  ^  >  14  km.  This  will,  at  the 
same  time,  enable  us  to  study  the  role  of  the  transition  region  in 
the  determination  of  the  transmission  and  reflection  coefficients. 
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In  Figs.  2  and  3.  we  have  plotted  It  I.  for  the  L-n«ode.  for 
srile  heights  of  ^  10  km  and  ^  »  l4  km,  respectively,  for  different 

values  of  the  ion  density.  (The  Gamma*- functions  of  complex  argument 
were  calculated  on  an  IBM  7044  electronic  computer.)  A  comparison 
of  these  figures  with  Fig.  1  for  the  layered  medium  shows  that  the 
most  striking  effect  of  the  inclusion  of  a  transition  region  is  the 
virtual  elimination  of  all  but  the  lowest  resonance,  which  is  in  turn 
considerably  depressed  and  broadened  compared  to  the  layered  case. 

This  qualitative  change  in  transmission  characteristics  results 
primarily  from  the  removal  of  the  artificial  sharp  boundary  of  the 
two-layer  limit.  The  lowest  resonance  still  remains,  because  at 
this  frequency  the  wavelength  is  many  time  the  scale  height  of  the 

I 

transition  region,  so  that  the  transition  from  Alfven  to  Hall  region 
appears  relatively  sharp.  At  the  higher  frequencies,  however,  the 
wavelength  is  con^arable  with  or  smaller  than  the  scale  height,  so 
that  one  region  appears  to  change  slowly  and  continuously  into  the 
other. 

The  effects  of  losses  in  the  transition  region  are  contained 

mainly  in  the  reflection  coefficient,  since,  with  our  model  of  a 

perfectly  reflecting  earth,  all  the  energy  which  is  not  reflected  is 

absorbed.  In  Figs.  4  and  5,  we  have  plotted  |r|,  for  the  L-mode, 

2 

for  three  different  densities.  (We  recall  that  |rj  is  the  energy 
reflection  coefficient.)  A  coi2q>arison  of  these  figures  with  Figs. 

2  and  3  shows  that  absorption  becomes  important,  for  a  given  ion 
density,  at  frequencies  somewhat  higher  than  the  corresponding  lowest 
resonance.  It  is  also  apparent  that  the  reflection  coefficient  ex¬ 
hibits  characteristic  minima  in  the  neighborhood  of  the  maxima  of  the 
corresponding  transmission  coefficient.  These  minima,  which  r'jpresent 
absorption  maxima,  become  more  pronounced  at  hi^er  frequencies. 

A  comparison  of  Figs.  2  and  3  with  Fig.  1  shows  that  the  intro¬ 
duction  of  a  transition  region  shifts  the  resonances,  but  not  drasti¬ 
cally,  to  somewhat  lower  frequencies.  These  resonances  are  thus 
essentially  still  those  frequencies  at  which  the  Hiill  region  becomes 
transparent.  The  downward  shift  produced  by  the  inclusion  of  a  transi¬ 
tion  region  results  from  the  fact  that  such  a  region  increases  the 
effective  size  cf  the  Hall  region. 
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Based  on  a  value  of  u)^  •  175  rad/sec,  the  lowest  resonances 
occur  at  a  frequency  of  .014  cps  for  ^  «  10  km  and  at  a  frequency 
of  .013  cps  for  2  -  14  km.  The  agreement  with  the  70  second  period 
observed  by  Ma^le  (1959)  is  undoubtedly  fortuitous,  but  the  general 
location  of  these  resonances  makes  It  seem  very  probable  that  they 
are  responsible  for  Pc  oscillations.  Another  point  In  support  of 
this  suggestion  Is  the  fact  that  these  resonances,  besides  being 
shifted  to  higher  frequencies  from  day  to  night,  are  also  considerably 
depressed  at  night.  This  would  account  for  the  fact  that  Pc  oscilla¬ 
tions  are  mainly  a  daytime  phenomenon. 

One  further  point  worthy  of  note  In  connection  with  Figs.  2  and 
3  Is  that  the  Ionosphere  becomes  essentially  transparen  t  night 
up  to  frequencies  of  the  order  of  1  cps.  On  the  other  .id,  in  the 
daytime  the  transmission  coefficient  becomes  rather  small  for  fre¬ 
quencies  greater  than  about  .05  cps.  This  would  account  for  the 
fact  that  "pearl"  t3rpe  oscillations,  which  have  frequencies  of  the 
order  of  .1  cps,  are  mostly  a  nighttime  phenomeno.i .  In  the  daytime, 
they  are  apparently  filtered  out  by  the  Hall  region. 

In  Figs.  6,  7,  and  8,  we  have  plotted  jt^j  for  the  R-mode,  l.e., 
for  (u  <  0.  It  Is  apparent  that  the  transmission  coefficient  for  this 
mode  drops  off  very  rapidly  with  frequency  compared  to  that  for  the 
L-mode.  At  resonance,  the  transmitted  wave  is  almost  entirely  of 
the  L-type,  and  consequently  left -circularly  polarized,  as  mentioned 
previously. 

In  Pigs.  9  and  10,  we  have  plotted  |r|  for  the  R-mode,  again 
for  three  different  ion  densities.  It  can  be  seen  that,  for  low 
frequencies,  this  mode  Is  almost  totally  reflected.  Indicating  that 
it  is  nearly  evanescent  at  these  low  frequencies.  For  frequencies 
higher  than  about  I  cps,  however,  significant  absorption  begins  to 
take  place,  especially  at  the  higher  Ion  densities. 

Finally,  we  should  mention  that  we  have  considered  only  the 
steady-state  problem.  The  transmitted  signal  resulting  from  a 
given  source,  however.  Is  determined  not  only  by  the  transmission 
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copf  f <  ri *»nf .  hut  also  bv  the  freouencv  snectruin  of  the  source.  For 

-  - - - - ^  ^  ^  ^  ^ 

a  given  source  spectrum,  however,  the  resulting  signal  can  be  calcula 
ted  from  the  transmission  coefficient  given  by  Eq.  (40).  This  in¬ 
volves  a  Fourier  integral  over  frequency  of  the  transmission  co¬ 
efficient  weighted  with  the  source  spectrum.  The  usual  technique 
for  evaluating  such  integrals  is  to  transform  them  to  contour  inte¬ 
grals  in  the  complex  ^-plane.  The  integral  then  becomes  the  sum  of 
the  residues  of  the  integrand  at  the  poles  of  the  transmission  co¬ 
efficient.  A  procedure  for  determining  the  poles  of  Eq.  (40)  in  the 
coiiq)lex  *u)-plane  has  been  presented  elsewhere  (Greif inger  and 
Gre if inger ,  1964).  At  this  point  we  would  only  remark  that  any 
transmitted  signal  becomes  the  sum  of  a  ntmuber  of  resonance-like 
terms,  with  coefficients  depending  on  the  source  spectrum.  The 
higher  frequencies  are  damped  out  more  rapidly  than  the  lower  ones, 
so  that  eventually  the  signal  is  almost  purely  sinusoidal,  its 
frequency  being  given  by  the  real  part  of  the  appropriate  pole  of 
the  transmission  coefficient.  This  is  very  nearly  the  value  of 
u)  at  which  the  transmission  coefficient  has  its  maximum  on  the  real 

uJ-exis,  i.e.,  at  the  resonant  frequency  calculated  above. 


Amplitude  of  transmitted  magnetic  field,  for  unit  amplitude 
of  incident  magnetic  field,  as  a  function  of  frequency,  for 
different  ionospheric  conditions.  (N.  is  the  ion  number 
density  in  units  of  10^/cc  and  a)^  is  the  ion  cyclotron 
frequency,) 


Anq  litude  of  transmitted  magnetic  field,  for  unit  anroli 
tuce  of  incident  magnetic  field,  as  a  function  of 
frequency,  for  different  ionospheric  conditions.  (Nj^ 
is  :be  ion  number  density  in  units  of  10^/cc  and  is 
the  ion  cyclotron  frequency.) 
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Amplitude  of  transmitted  .iUgnetic  field,  for  unit 
amplitude  of  incident  •r.agnetic  field,  as  a  function 
of  frequency,  for  different  ionospheric  conditions. 
(Ni  is  the  ion  number  density  in  units  of  10^/cc  and 
U),-  is  the  ion  cyclotron  frequency.) 
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Amplitude  of  reflected  magnetic  (or  electric)  field,  for 
unit  amplitude  of  incident  magnetic  (or  electric)  field, 
as  a  function  of  frequency,  for  different  ionospheric.  ^ 
conditions.  (N^  is  the  ion  number  density  in  units  of  lO”/cc 
and  is  the  Ion  cyclotron  frequency  ) 


Aniplitude  of  reflected  magnetic  (or  electric)  field,  for 
unit  amplitude  of  incident  magnetic  (or  electric)  field, 
as  a  function  of  frequency,  for  different  ionospheric 
conditions.  is  the  ion  number  density  in  units  of 

10^/cc  and  ai.  is  the  ion  cyclotron  frequency.) 


o.o<^ 

\ 

\ 

km 

(R-Mode) 


Fig.  6  -  Amplitude  of  transmitted  magnetic  field,  for  unit 

amplitude  of  incident  magnetic  field,  as  a  function 
of  frequency,  for  different  ionospheric  conditions. 
(Nji  is  the  ion  number  density  in  units  of  10^/cc  and 
u).  is  the  ion  cyclotron  frequency.) 
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Fig.  7  -  Amplitude  of  transmitted  magnetic  field,  for  unit 

amplitude  of  incident  magnetic  field,  as  a  function  of 
frequency,  for  different  ionospheric  conditions.  (N^ 
iS  the  ion  number  density  in  units  of  10^/cc  and 
is  the  ion  cyclotron  frequency.) 
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Amplitude  of  transmitted  magnetic  field,  for  unit 
amplitude  of  incident  magnetic  field,  as  a  function 
of  frequency,  for  different  ionospheric  conditions. 
(Ni  is  the  ion  nuir.ber  density  in  units  of  10^/cc  and 
uu.  is  the  ion  cyclotron  frequency.) 


Fig.  9  -  Anq>litude  of  reflected  magnetic  (or  electric)  field, for 
unit  amplitude  of  incident  magnetic  (or  electric)  field 
as  a  function  of  frequency,  for  different  ionospheric 
conditions.  (N^  is  the  ion  number  density  in  units  of 
10^/cc  and  uu^  is  the  ion  cyclotron  frequency.) 


Amplitude  of  reflected  magnetic  (or  electric)  field, 
for  unit  amplitude  of  incident  magnetic  (or  electric) 
field,  as  a  function  of  frequency,  for  different 
ionospheric  conditions.  (N.  is  the  ion  number  density 
in  units  of  10°/cc  and  x,  is  the  ion  cyclotron  frequency.) 
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APPFMnTY 


In  order  to  solve  Eq .  (24),  we  make  the  following  change  of 


variable : 


1  ,  -z/r  ,  ^  -i^/2  -zfl 

u=-l-ie  ^  =  l-»-e  e 


In  teims  of  this  variable,  Eq.  (24)  becomes 


2^2 


t-l)du  u(u-U  J 


E  =  0, 


This  is  to  be  compared  with  the  canonical  form  of  the  hyper geometric 
equation 


2  L  u  (u-1)  J  du  L  ^"l)  j  u(u-l) 


These  equations  become  identical  if  we  set 
a  *  -o'  •  i^  (u 

6  -  -  0 '  -  i'J'  u)  ^  (A4 

*y  -  0,  y'  •  1  • 

Ue  are  looking  for  solutions  of  Eq.  (A2)  which  have  the  form 
prescribed  by  Eq.  (28),  viz.. 


-  e 


W  -•  e 
“2  ® 


itt)  z 


-ito  z 


(z  -•  ») . 


In  terms  of  the  variable  u  defined  by  Eq.  (Al) ,  these  solutions 
must  have  the  form 
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„  n/2  Z  ^  ^  e 

Wj^  -  e  (u-1) 


(u  -  1)  . 


(A6) 


The  solutions  of  Eq.  (A2)  which  have  the  form  required  by  Eq.  (A6) 
are 


TT  ^  ^ 


w. 


UU  "  ,-v  ^ 


“  F(-i?  S  +  1?  I.  -  If  S  ^  S; 


w. 


1  -  2ir  uJ  -  (U-D) 

e  (u-l)^^*"  F{iZlS^  -t-  iZH,  iZ'^^  -  iZ'^i 


(A7) 


1  4-  2i'c  w  -  (u-1)) , 


where  F  is  the  usual  hyper^ceonetric  function. 

We  now  need  to  evaluate  Wj^  and  W2  in  the  limit  z  which, 

Trom  Eq.  (Al)  ,  is  the  limit  u  -«  ob.  To  do  this,  we  make  use  of  the 
analytic  continuation  which  relates  the  hypergeometric  function  of 
arguoient  x  and  the  hyper  geometric  function  of  argument  x  viz.. 


F(a,b;c;x)  -  ^r(c-a)  1-b  +  a;  x  ^ 


+  C-*)'**  F(b,  l-c4-b;  1-a+b;  x'b 

I  arg  (-X)  I  <  TT, 


(A8) 


where  T  is  the  usual  r^function.  We  apply  Eq.  (A8)  to  Eq.  (A7)  and 
take  the  limit  u  -*  es,  where  the  hyper geosK trie  functions  may  be 
replaced  by  unity.  It  is  then  easy  to  show  that  Eq.  (A7)  becoaies 


1 

2 


iu)  Z  .  .  -Itt)  z 

a^  e  +  bj^  • 

iUB  z  .  .  -i»  z 

^2  +  *>2  * 


(z  -  -  •) 


(A9) 
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i^  'z 

on  the  right-hand  side  of  Eq .  (A7) ,  while  the  outgoing  wave,  e  , 
arises  entirely  from  the  second  term.  The  coefficients  in  Eq . 

(A9)  are  given  by 

2  ^  ^  ^  r(i  -  21?  ^  S  r(-2i^;a _ 

r(-i?  ^  i'C  ^  r(i-i^:  ^  ^  -  i'^:  ^ 


ai 


.  e 


-  ?r  ® 


-  ?  r  3  ^  r(i  >  ai?  s  S  r(-2ir  3} 


-  ?rs 


a„  »  e 


r(ic  00 


ir  3)  rcn-  IT  5  ^  -  1?  3) 


(AlO) 


rti-atrs^  r(2irs> 


•?t 


U) 


r(-i^  iZ'SS)  rci-i'f  ^  ^  +  ir  ^ 


.  e 


■  2^“  ■  r(i  21? s r(2tr a _ 2 

r(i73;  *  +  i7S)  r(i+  i^fs  V  i? s) 


“rs 


To  obtain  Eqs.  (40)  and  (41),  we  substitute  Eq.  (AlO)  into 
Eqs.  (36)  and  (38),  respectively,  and  make  use  of  the  recursion 
formula  r(z+l)  *  z  r(z) . 
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